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Abstract We analyze the optical spectrum of type 1 QSO SDSS J1425+3231. This ob- 
ject is interesting since its narrow emission lines such as [O III] AA4959, 5007 are double- 
peaked, and the line structure can be modeled well by three Gaussian components: two 
components for the two peaks (we refer the peaks at low/high redshift as "the blue/red 
component") and another one for the line wing which has the same line center as that 
of the blue component, but ~ 3 times broader. The separation between the blue and red 
components is ~ 500 km/s with blue component ~ 2 times broader than the red one. The 
H/3 emission can be separated into four components: two for the double-peaked narrow 
line and two for the broad line which comes from the broad line region (BLRs). The black 
hole mass estimated from the broad H/3 emission line using the typical reverberation map- 
ping relation is 0.85 x 10 8 M Q , which is consistent with that derived from parameters of 
[O III] A 5007 of the blue component. We suggest this QSO might be a dual AGN system, 
the broad H/3 emission line is mainly contributed by the primary black hole (traced by 
the blue component) while the broad H/3 component of the secondary black hole (traced 
by the red component) is hard to be separated out considering a resolution of ^2000 of 
SDSS spectra or it is totally obscured by the dusty torus. 
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1 INTRODUCTION 

Early in 1980s, the double-peaked narrow emission lines in active galactic nuclei (AGNs) have been 
reported by Heckman et al. (1981, 1984), and they suggested biconic outflows as an origin of double 
peaked AGNs. Greene & Ho (2005) pointed out that there are about 1% local AGNs have double-peaked 
narrow emission lines. Since narrow emission lines are generally believed to be produced by clouds in 
the narrow line regions (NLRs), they suggested that ^1% AGNs have disk-like NLRs as a simplest 
explanation. The ~1% double-peaked AGN fraction have also been found by Zhou et al. (2004) and 
Wang et al. (2009). Especially, the flux ratios of the two peaks of [O III] A 5007 are anti-correlated 
with the ratios between their shifts relative to the host galaxies statistically, this phenomenon could be 
understood as a natural result of Kepler rotation of dual AGNs, pointing to another possible explanation 
of double-peaked narrow emission lines (Wang et al. 2009). 

In the last two years, more and more attention has been paid to the study of dual AGNs (e.g. 
Comerford et al. 2009a, 2009b; Dotti et al. 2009; Boroson & Lauer 2009; Xu & Komossa 2009; Liu 
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et al. 2010) since on the one hand, they are unavoidable results of hierarchical cosmology model which 
has gained great success; on the other hand, with the launch of X-ray and infrared space astronomical 
facilities (Chandra, Spitzer, and Herschel Space telescopes), it is now possible for us to resolve the dual 
AGNs with kpc-scale separation. Observational evidence for dual AGNs includes: spatially resolved 
systems in which both supermassive black holes (SMBHs) can be identified directly and spatially unre- 
solved systems in which the dual AGN model can explain various phenomena (see Komossa 2006 for 
a detail review). Most recently, Colpi & Dotti 2009 further summarize the observations and numerical 
simulations of dual and binary black holes. So far, a few unambiguous cases have been found, such as 
NGC 6240 (Komossa et al. 2003), J0402+379 (Rodriguez et al. 2006, 2009; Morganti et al. 2009), 
EGSD2 J142033+525917 (Gerke et al. 2007), EGSD2 J141550+520929 (Comerford et al. 2009a), 
COSMOS J 100043+020637 (Comerford et al. 2009b), and other four dual AGNs (see Liu et al. 2010). 
There are other interesting sources, which need more certification by future observations. Zhou et al. 
(2004) connected both SDSS and VLBA data, suggesting SDSS J1048+0055 is a dual AGN system and 
double-peaked narrow emission lines could be an effective way of selecting dual AGN candidates. Xu 
& Komossa (2009) analyze the line structures and flux ratios of SDSS J13 16+1753 in detail, discussing 
all the possible origins of the double peaks. Furthermore, SDSS J1536+0441 is the only source in which 
two broad line systems have been found, and it is suggested to be a binary black hole system which 
is separated by O.lpc with a orbital period of 100 yr (Boroson & Lauer 2009). However it is also the 
most controversial case, which brings about a vast deal of debate (Chornock et al. 2009, 2010; Wrobel 
& laor 2009; Decarli et al. 2009; Lauer & Boroson 2009; Tang & Grindlay 2009; Gaskell 2010; Dotti & 
Ruszkowski 2010; Bondi & Perez-Torres 2010). 

The peculiar emission-line spectrum of SDSS J1425+3231 is noticed in the course of searching for 
dual AGN candidates in the SDSS QSO sample. SDSS J 1425+3231 shows all its strong narrow emission 
lines with double-peaked, it is hard to make conclusions on the structures of weak emission lines due to 
the signal-to-noise ratio (S/N). We analyze the spectra of SDSS J 1425+3231 in §2. The black hole mass 
is estimated in §3. The possible origins of the double-peaked line profiles of this source are discussed in 
§4. The results are summarized in §5. Throughout this paper, a cosmology with Hq= 70 km s _1 Mpc _1 , 
f2jvf= 0.3, and J7a= 0.7 is adopted. 

2 DATA ANALYSIS 

SDSS J 1425+3231 is a broad line QSO with SDSS pipeline redshift of z = 0.478, all its strong narrow 
emission lines show double-peaked profiles. In this section, we describe the procedure of spectral fitting. 
The steps of our analysis are as follows: (1) the spectrum is corrected for foreground Galactic extinction 
and shifted to the rest-frame by using z = 0.478; (2) the continuum of the spectrum is modeled by three 
components (Hu et al. 2008) and subtracted, the aim is to separate the contribution of continuum and 
emission line spectrum; (3) multiple Gaussian components are used to fit the emission lines. Step (2) 
and (3) are described in more detail below. 

2.1 Continuum decomposition 

The continuum is modeled as 

F\ = F A PL (F 510 o, a) + F A BaC (F BE , r BE ) + F* c (F Fc , FWHM Fo , V Fc ). (1) 

where — -FiooCi^o)" is a featureless power law, F5100 is the flux at 5100 A and a is the spectral 
index. The second and third terms represent the Balmer continuum and Fe emission, respectively. 

For wavelength shortward of the Balmer edge A < Abe = 3646A, the Balmer continuum can be 
expressed as Ff aC = F B eB x (T e )(l- e~ TA ) (Grandi 1982; Dietrich et al. 2002). F B e is a normalization 
coefficient for the flux at Abe, Bx(T e ) is the Planck function at an electron temperature T e , tx = 
t~be(j^;) is the optical depth at wavelength A, tbe is the optical depth at the Balmer edge. T e is 
assumed to be T e = 15,000K. The two free parameters in the Balmer continuum are Fbe and tbe- 
At A > Abe, blended higher-order Balmer lines give a smooth rise from ^4000 A to the Balmer edge 



SDSS J142507. 32+323 137.4 



3 



(Wills et al. 1985) in the spectrum. However, our fitting windows do not include this region, actually our 
results are not influenced by the higher order Balmer lines. 

The optical and ultraviolet Fe II template (Fj[ Zwl ) from NLS1 I ZW 1 is used to subtract the Fe II 
emission from the spectra (Boroson & Green 1992; Vestergaard & Wilkes 2001). The I ZW 1 template 
is broadened by convolving with a Gaussian function G: 

F Fe = F IZwl „ G (F conv ,FWHM conv , Konv), (2) 

where F conv , FWHM conv , and V^ onv are the flux, width and peak velocity shift of the Gaussian function. 
The parameters of Fe in equation (1) can be expressed as follows: the flux of the Fe emission, Fp e , it is 
the multiplication of F conv and the flux of the template; the shift of the Fe spectrum, Vp c = V^ onv , and 
the FWHM of the Fe lines FWHM Fc = v /FWHMj ! Zwl + FWHM^ nv . 

In total, there are seven parameters in the continuum model, they are fitted by minimizing x 2 - 
The fitting windows include: 2470-2625, 2675-2755, 2855-3010, 3625-3645, 4170-4260, 4430-4650, 
5080-5550, and 6050-6200 A. These windows are free of strong contaminant lines. Figure 1 shows 
the result of the continuum decomposition. The Galactic extinction and redshift corrected spectrum is 
shown in the top panel. The spectrum in the fitting window is plotted in green. The three components 
of the continuum are shown in blue. The best fit model is shown in red. The middle panel shows the 
residual spectrum, namely the pure emission-line spectrum. We will analyze it in the next step. In the 
bottom panel, we subtract the power law and the Balmer continuum, zoom in Fe-only spectrum in the 
wavelength range 4100-5600 A. The Fe model is shown in red. 

2.2 Emission line fitting 

We measured the H/3 and [O III] A4959, 5007 emission lines from the emission line spectrum. Since 
the narrow lines of this object are double-peaked, two sets of three Gaussian components are used to 
model [O III] A4959 and 5007. One for the blue component, one for the red component, and one extra 
Gaussian component for the underlying broad wing. Fit parameters of all the three Gaussians are the 
central wavelength, a and flux. Each [O III] A4959 component is forced to have the same centroid and 
width as the corresponding [O III] A5007. In Figure 2, the blue and red components are shown in blue 
and red respectively. The orange is the underlying wing. The best fit model is shown in green. The 
three-Gaussian parameterization describes the observed [O III] A4959, 5007 lines very well. The blue 
and red components are separated by ~500kms _1 in velocity space, and the a of blue component is 
417km s~\ about 2.2 times broader than the red component. The broad wing (orange) peaks at roughly 
the same position as that of the blue component. The a of the wing is 628km s~\ We use four-Gaussian 
components to model the asymmetric profile of the H/3 emission, two for the blue and red narrow peaks, 
they are forced to have the same width as the corresponding [O III] A5007, see the blue and red Gauss 
in Figure 2. Different from the fit of [O III], we do not fit a wing to the narrow H/3 component since 
on the one hand, comparing with the broad H/3, the strength of the wing can be neglected; on the other 
hand, the resolution of the SDSS spectrum is not high enough to separate such a weak component from 
the broad H/3. The broad H/3 are fitted with the two other Gaussian components, see the pink lines. The 
fitting parameters of each Gaussian component is shown in table 1 . 

3 MASS OF THE CENTRAL BLACK HOLE 

In this section, we use three different methods to estimate the mass of the central black hole. 

The first method we used is based on the virial theorem and -Rblr — £5100 relation which is cal- 
ibrated with the reverberation mapping data. The BLR radius -Rblr is estimated from the continuum 
luminosity at 5100A (£5100 = 1-6 x 10 40 erg s _1 A -1 ) using the -Rblr — £5100 relation given by 
Bentz et al. (2006), then the black hole mass is estimated from A/bh = f Ssi-^Z. ; where / is the 
scaling factor which is introduced to characterize the unclear kinematics and geometries of the BLRs. 
Actually, the value of / changes with the shape of the line in use, / = 3.85 is a mean value suggested 
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Fig. 1 SDSS spectrum of SDSS J 1425+3231 plotted as observed flux versus rest wavelength. 
The top panel shows the Galactic extinction and redshift corrected spectrum. The spectrum in 
the fitting window is plotted in green. The three components of the continuum are shown in 
blue while The best fit model is shown in red. The middle panel shows the residual spectrum, 
namely the pure emission-line spectrum. The bottom panel shows the spectrum after subtracting 
the power law and the Balmer continuum in the wavelength range 4100-5600 A. The Fe model 
is shown in red. The flux density F\ is given in units of 10~ 17 erg s _1 cm _2 A _1 



by Collin et al. (2006). AV — 1325km s 1 is the second moment of the BLR H(3 profile, which is 
reconstructed from the two pink components in Figure 2. This method gives a value of 8.5 x 10 7 A'/ o . 

The Mbh — £5100 relation (see Eq. 9 of Peterson et al. 2004) predicts a mass of 3.50 x lO 8 A/0. 
We also use velocity dispersion of gas in the NLR, namely a of [O III] A 5007 as a surrogate of galaxy 
stellar velocity dispersion c* (Nelson 2000), to estimate the central black hole mass based on the famous 
M Ba - o* relation which is in the form of log(M B H/M ) = 8.13 + 4.021og(cr„/200)(Tremaine et al. 
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Fig. 2 Emission-line fitting of SDSS J1425+3231. For [O III] A4959, 5007 emission lines, three 
Gaussian components are used. The blue and red Gauss represent the two narrow line peaks, the 
orange is the underlying wing. The broad Hf3 are fitted with four Gaussian components, two for 
the double-peaked narrow line (red and blue) and another two for the broad line (see the two pink 
components). The best fit model is shown in green. 



2002). The blue component gives a mass of 8.33 x 10 7 A/q with a = 177.40 km s 1 while the red 
component indicates a mass of 3.42 x 10 6 Af Q with a = 80.18 km s _1 . 



4 DISCUSSION 

4.1 Superposition of two objects 

This possibility is very unlikely for two reasons. First, the red component of [O III] A 5007 has a luminos- 
ity of 3.2 x 10 42 erg s _1 , while the a of this line indicates a black hole mass of 3.42 x 10 6 M©. If this 
component comes from a background AGN, its accretion rate should be extremely super-Eddington, 
which is hard to explain (see section 4.3 for the explanation under the dual AGN scenario). Second, 
Dotti & Ruszkowski (2010) examined the superposition model of double-peaked emission line AGNs 
based on galaxy clusters from the Millennium Run, finding that the fraction of superimposed galaxy 
pairs peaks at about z = 0.2 and decreases rapidly since z = 0.3. Considering the redshift of SDSS 
J1425+3231, ~ 0.478, the possibility of superposition should be very low. 
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Table 1 Emission line properties of SDSS J1425+323 1 



Property H/3 [O III] [O III] 
4861.33 4958.91 5006.84 

Blue system 

Line center (A) 4861.77±0.51 
FWHM (km s _1 ) 417.17±17.77 
Luminosity" 44.51±9.88 

Red system 
Line center (A) 
FWHM (km s _1 ) 
Luminosity 11 

Wing 

Line center (A) 
FWHM (km s" 1 ) 
Luminosity" 

Broad system 
Line center (A) 
FWHM (km s" 1 ) 
Luminosity" 
Very broad system 6 
Line center (A) 
FWHM (km s" 1 ) 
Luminosity" 

Notes: 

a In units of 10 40 erg s — 1 . 
b Note only the H/3 need the very broad component in appearance. 

4.2 NLR kinematics 

Another possible explanation for the double-peaked narrow emission could be special NLR geometries 
such as biconical outflows and disk-like NLRs. In such hypothesis, there is only one AGN to illuminate 
the NLR gas which is moving toward and away from us, forming the blue and red components in the 
observed narrow lines. 

Certain nearby Seyfert and star forming galaxies are known to have biconic outflow induced double- 
peaked emission lines. The examples are found not only from the spatial resolved spectra which takes 
along the minor axis of a galaxy (e.g. Cecil 1988; Cecil et al. 1990; Veilleux et al. 1994, 2001; Colbert 
et al. 1996) but also from the spectrum of the whole galaxy (e.g. Duric & Seaquist 1988; Axon et al. 
1998). In the scenario of biconic outflows, we would expect that the blue and red components have the 
same velocity dispersion since they are illuminated by the same AGN. However, in SDSS J 1425+3231, 
the blue component is three times broader than the red component, which is conflicted with the outflow 
model. On the other hand, the outflow studies find that the NLRs are stratified strongly in ionization and 
velocity so that high-ionization lines, such as [O III] are originated near the AGN with higher velocity 
and low-ionization lines, such as H/3 and [O II], are originated further from the AGN with lower velocity 
(e.g., Komossa et al. 2008). In SDSS J1425+323 1, we find that [O III] and Hf3 are consistency in velocity 
offsets within error bars, and apparently no ionization stratification is observed as expected in AGN 
driven outflows. In addition, the [O lll]A5007 luminosity is 4.3 x 10 42 erg s _1 for the blue system and 
2.8 x 10 42 erg s _1 for the red system. These are typical values for the emission from the whole NLRs 
of bright Seyfert galaxies and quasars, but it is hard to image that the whole NLRs are outflowing. 

The observed double-peaked narrow emission lines can be accounted for in a disk-like NLR model. 
The rotating disk model predicts that the blue and red components have similar width, which is conflict 
with the current data. Furthermore, the red and blue components are expected to be (almost) equally 
shifted with respect to the true cosmological galaxy redshift in this scenario. We are lack of host galaxy 
information since it is over-shined by the central AGN, however under the disk-like NLR model, the 



4957.59±0.08 
417.17±17.77 
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broad line should be at the redshift of the galaxy, in between the blue and red components, and not so 
close to the the blue component as we have seen from the data. 

4.3 A dual AGNs system 

The final picture we want to suggest for SDSS J 1425+3231 is that this is a dual AGN system. The 
primary black hole in this system has a mass of ~ 10 8 Af Q , the observed blue narrow component in 
[O III]AA4959, 5007, the broad H/3 component and the wing of [O lll]AA4959, 5007 are generated in 
its NLR, BLR and the region in between (this region is usually referred as intermediate line region 
in literatures) respectively. The consistency between the black hole masses estimated from broad H/3 
component and a of the blue component supports this idea. We have not observed the broad emission 
lines from the BLR of the secondary black hole (~ 10 Mq), the possible reasons for this includes: (1) 
the S/N of SDSS spectra is not high enough for us to separate the broad emission line of the secondary 
black hole from primary one; (2) the secondary black hole is a type 2 AGN in which the BLR is obscured 
by the dusty torus. 

The only issue about the dual AGN system is that, for the secondary black hole, comparing with the 
[O III] A 5007 luminosity, a mass of 10 6 Af© is a little bit low. This indicates the secondary black hole is 
accreting in a super Eddington regime. At first sight, this is conflicted with the normal accretion theory. 
However, we should note that in the dual AGN case, the separation between the two black holes is in 
kpc scale, the NLR gas of the secondary black hole can be affected by the primary black hole, namely 
the primary black hole can illuminate the NLR of the secondary black hole. 

5 CONCLUSION 

Type 1 QSO SDSS J1425+3231 has double-peaked narrow emission lines. In this paper, we analyze 
the SDSS spectrum of this object, discussing the origins of its double-peaked line structure. We argue 
against the possibility of superposition of two objects, biconic outflow and disk like NLR, suggesting 
that this is a dual AGN system. 

In this system, the primary black hole has a mass of ~ 1O 8 M , the observed broad lines belong to 
it. The secondary black hole is much smaller, its mass is in the oder of 10 6 M Q . The secondary black 
hole could be a type 2 AGN whose BLR is obscured, or we failed to separated the broad emission line of 
the secondary black hole from that of the primary black hole due to the resolution of the SDSS spectra. 

In the future observation, high spatial resolution two-dimensional optical spectrum, and imaging 
in the optical, radio, X-ray would help us figure out whether SDSS J 1425+3231 contains dual AGNs. 
Moreover, the ultimate confirmation or rejection of the dual AGN interpretation, which predict the vari- 
ations in the line profiles, will come from multi-epoch spectroscopic monitoring. 
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